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Abstract 



^> . We demonstrate that parastatistics can be quantized using path integrals 

\D ■ 

by calculating the generating functionals for time-ordered products of both 



free and interacting parabose and parafermi fields in terms of path integrals. 



We also give a convenient form of the commutation relations for the Green 



components of the parabose and parafermi operators in both the canonical 
and path integral formalisms. 



1 Introduction 

Canonical quantization [1, 2] and path integral quantization [3, 4] are two widely 
used methods to quantize a physical system. The generality of these two methods 
and their equivalence or non-equivalence has general interest. Many physicists, 
including the present authors, believe that any system that can be quantized using 
one of these methods also can be quantized using the other. We acknowledge that 
for technical reasons one or the other method may be simpler or more convenient 
for any specific system. 

Some years ago M.G.G. Laidlaw and C. M. DeWitt [5] in a well-cited influential 
paper studied path integral quantization for a system of identical particles. They 
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noted the important distinction between quantization of identical particles in two 
space dimensions and in three or more dimensions. This distinction holds for both 
path integral and canonical quantization. We believe that their paper correctly 
proves that parastatistics cannot be based on the topology of three dimensional 
space (in contrast to anyon statistics [6, 7], which is based on the topology of two 
dimensional space). Although the text of their paper does not mention parastatis- 
tics, their abstract contains the phrase "... showing that the Feynman formalism 
excludes parastatistics." Dr. Cecille DeWitt has informed us that the abstract 
is misleading and that M. Laidlaw and C. DeWitt did not intend to state that 
there is a contradiction between path integration and parastatistics [8]. If, indeed, 
parastatistics fields [9] could not be quantized using path integrals that would have 
significance for the general question of the equivalence of the canonical and path 
integral methods of quantization. Our belief that any field theory that can be quan- 
tized using the canonical formalism also can be quantized using path integrals was 
our initial motivation to formulate field theories with parastatistics fields using path 
integrals. 

There are several plausible assumptions about the quantum mechanics of in- 
distiguishable particles that lead directly to the conclusion that only bosons or 
fermions can exist in three or more space dimensions. These plausible assumptions 
lead to the same conclusion in the context of path integrals. The basic principles 
of quantum mechanics do not require these plausible assumptions; indeed without 
these assumptions quantum mechanics allows the description of indistinguishable 
particles that have more general statistics than bose or fermi statistics [10, 11]. 
Here are some of the plausible assumptions together with the reasons why these 
plausible assumptions are not required for quantum mechanics: (1) Transposing 
indistinguishable particles in a quantum state produces a state that is linearly de- 
pendent on the initial state, i.e. T\& = rjty, where T is a transposition. Since the 
square of a transposition is the identity this leads directly to T 2 \l/ = rf^> = \&, so 
r\ — ±1. In a general representation of the symmetric group this is not the case. 
Greenberg and Messiah showed [10] that quantum mechanics can be formulated for 
an arbitrary represention of the symmetric group. (2) The space of quantum states 
obtained by permuting indistiguishable particles is coherent, i.e., these states can 
be superposed to form other states without restriction. For states of several indis- 



tinguishable particles there is an absolute selection rule that prohibits transitions 
between states in inequivalent irreducible representations of the symmetric group 
so that the relative phases of states in inequivalent irreducible representations can 
never be measured. Such systems must be treated as statistical mixtures rather 
than as coherent superpositions of states. In short, there is a superselection rule 
that separates states in inequivalent irreducible representations of the symmetric 
group [10, 12]. If one assumes all states can be superposed, then since there is a 
symmetric state and an antisymmetric state of N particles for each N, only bosons 
and fermions can occur. (3) The configuration space of a system of particles is in 
one-to-one correspondence with the states of the system. This assumption implies 
that the states always can be superposed and leads back to assumption (2). (4) 
The homotopy classes of the configuration space determine the particle statistics. 
The discussion of irreducible representations of the symmetric group shows that the 
existence of different homotopy classes allows different types of statistics, but there 
is no requirement that each type of statistics must be associated with a different 
homotopy class of the configuration space. 

2 Canonical Quantization of Parastatistics 

H.S. Green [9] proposed the first proper quantum statistical generalization of bose 
and fermi statistics. Green noticed that the commutator of the number operator with 
the annihilation and creation operators is the same for both bosons and fermions 

[n k ,aj]- = Skia}. (1) 

The number operator can be written 

n k = (l/2)[4,a fc ]± + const., (2) 

where the anticommutator [aj.,Ofc]+ (commutator [a[,Ofe]_) is for the bose (fermi) 
case. If these expressions are inserted in the number operator-creation operator 
commutation relation, the resulting relation is trilinear in the annihilation and cre- 
ation operators. Polarizing the number operator to get the transition operator n k i 
which annihilates a free particle in state I and creates one in state k leads to Green's 



trilinear commutation relation for his parabose and parafermi statistics, 

[[4> a j]±»4»]- = 2 <*Jm4 (3) 

Since these rules are trilinear, the usual vacuum condition, 

a k \0) = 0, (4) 

does not suffice to allow calculation of matrix elements of the a's and a^'s; a condition 
on single-particle states must be added, 

a k a\\0)= P 5 kl \0). (5) 

Green found an infinite set of solutions of his commutation rules, one for each 
positive integer p, by giving an ansatz which he expressed in terms of bose and fermi 
operators. Let the annihilation and creation operators be represented by a sum of 
operators with an additional index, the "Green" index, 



a 



l = tt )] ^'■ l '° 



b£ n , « t = Ef, (6) 



a=l q=1 



and let the b^ and bf be bose (fermi) operators for a = (3 but anticommute 
(commute) for a ^ (3 for the "parabose" ("parafermi") cases. This ansatz clearly 
satisfies Green's relation since the (polarized) number operator is diagonal in the 
Green index for each case. The integer p is the order of the parastatistics. The 
physical interpretation of p is that, for parabosons, p is the maximum number of 
particles that can occupy an antisymmetric state, while for parafermions, p is the 
maximum number of particles that can occupy a symmetric state (in particular, the 
maximum number which can occupy the same state). The case p — 1 corresponds 
to the usual bose or fermi statistics. Later, Greenberg and Messiah [10] proved that 
Green's ansatz gives all Fock-like solutions of Green's commutation rules. 

We give a useful form of the commutation relations for the Green components 
of a parabose (parafermi) operator, 



f )f = ±(25 Qi/3 - l)bi m b { k a) + SajSv, (7) 



where the upper (lower) sign is for parabose (parafermi); analogous forms hold when 
both operators are annihilation or both are creation operators, but the 8 a ,p$k,i term 
is absent. As far as we know this form does not appear in the literature. We find 
that the ±(2S a ^ — 1) factors are very useful in doing systematic calculations. 

We discussed quantization of parastatistics above in terms of the commutation 
relations for parabose and parafermi creation and annihilation operators. Quanti- 
zation of parastatistics in terms of spacetime-dependent field operators follows for 
free fields or fields in the interaction picture by the usual Fourier transformations 
using, for example, plane waves. Thus for a neutral scalar parabose field, 

m = j^j- 2 J -^Ifoe-** + ay k -% u k = V^Tml. (8) 

In terms of fields, local observables have a form analogous to the usual ones; for 
example, the local current for a spin-1/2 parafermi theory is 

U = (l/2)[^(x), 7 ^(x)]_ = J2^ a \x) 1 ^ a \x). (9) 

a=l 

From Green's ansatz, it is clear that the squares of all norms of states are positive, 
since sums of bose or fermi operators give positive norms. Thus parastatistics gives 
a set of orthodox theories. 

We give the form of the field commutation relations for parabose and parafermi 
fields in terms of both the the fields and the Green components of the fields, 

Mx 1 ), 0(x 2 )]±, 0(x 3 )]_ = 2iA(x 2 - x 3 )0(xi) ± 2(f ) (x 2 )iA(x 1 - x 3 ) (10) 

ct>^\x l )^\x 2 ) = (26 a ,f, - l)<P^(x 2 )^ a \x l ) + 6 a ,piA( Xl - x 2 ) (11) 

For charged fields, keep only the contractions between a field and its adjoint. 

3 Green Ansatz for Path Integrals 

3.1 Parabose case 

Since parafields can be constructed in second quantized form using the Green ansatz [9], 
we have constructed the representations in terms of path integrals by using construc- 
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tions closely analogous to those of Green. Since, as mentioned above, Greenberg and 
Messiah [13] proved that any field theory that obeys the trilinear commutation rela- 
tion of Green that has a Fock-like representation can be represented using the Green 
ansatz, the Green ansatz for the path integral can be used without loss of generality. 

We recall the generating a (bose) charged scalar field given in terms of a path 
integral, 



w[j, r 



I V<p* V<f> exp 



i / d 4 x (C((f)(x), <j>*(x)) + J*(x)(f)(x) + <j>*(x)J(x)) 



(12) 



with 

C = C + C t , C = d^d^ - m 2 0>, d = V(<f>, 0*). (13) 

Since one can calculate the path integral for the free term Co explicitly, one often 
represents the interaction term Ci by a functional derivative. Then 

w[j,r] = »'Wh[i/^(-i^,-^)l 

x exp[-i f d 4 xd A y.r(x)A F (x - y)J(y)], (14) 



where Ap(x) = (2tt) 4 / d 4 k exp(—ik ■ x) (k 2 — m 2 — ie) 



-i 



We construct the generating functional for time ordered products for parabose 
fields in parallel with the above construction for bose fields. We must take into 
account the requirement that the terms allowed in the parabose (and also in the 
parafermi) Lagrangian must obey local commutativity [13, 14]. These papers show 
that the para fields must occur in one of two forms. The form [<j)*(x), 4>(x)]± for any 
order p of parastatistics is analogous to the requirement that charged fields occur 
in the combination (f>*(x)(f>(x) in order to obey a U(l) symmetry. The form of a 
nested set of p — 1 commutators and anticommutators for para fields of order p is 
analogous to the possibility of a term e a,l3 '"'^(f> a {x)(j>f3(x) ■ ■ ■ 4>((x) in a theory that 
obeys SU(n) symmetry. We refer to the papers just cited for a demonstration and 
detailed analysis of these possibilities. In the first case the admissible terms in the 
Lagrangian will be functions of terms diagonal in the Green index. In the second 
case the admissible terms will be functions of terms antisymmetric in the Green 
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index; that is for a given value of p each index from 1 to p will occur exactly once 
in each term. 

The Lagrangian density for a free parabose field is 

£(0,0*) = \([d^(x),dy(x)} + -m 2 [0*(a;),0(x)] + ). (15) 

Using the Green ansatz, the parabose field is represented by 

#s) = X> (a) (s), (16) 

a=l 

that is, as a sum of p bose fields (p is the order of the parastatistics). The bose fields 
with different Green index anticommute rather than commute. In order to obtain 
the generating functional, we introduce parabose Grassmann fields J(x) [15]. The 
parabose Grassmann fields have a representation similar to that for the parabose 
fields, except that the right-hand sides always vanish, 



J(x) = Y,J {a) (x), (17) 

a=l 

[J^(x),J^(y)} qaf3 = 0, (18) 

where q a p = 1 — 25 a p. The parabose Grassmann Green components commute with 
the parabose field Green components for the same value of the Green index and 
anticommute for different values of the Green index. The external source term 
J d 4 x(l/2)([J*(x), <f>(x)] + + [J(x),<f>*(x)]+) is then an effective bose operator [15]. 
Using the fact that 

hd^(x),dy(x)] + = -£d^*(x)d^ a \x), (19) 

1 i 

we write the generating functional for a free parabose field as 

f[ V<p {a) *V<p {a) exp[-i^ fd 4 x(d^ {a) *(x)d'"(f) {a) (x) -m 2 <p {a) \x)<P {a \x) 

a=l 1 

+ {J {a) \x)^ a \x) + h.c.))}. (20) 



As in the bose case we can do the free path integral to find 

W free [J,J*] = W free [0,0]exp[-iJ2jd 4 xd 4 yjW*(x)A F (x-y)j( a \y)]. (21) 



a=l' 



For interacting theories in which the terms are functions of [<f)*(x),4>(x)]±, 
the interaction term in the Lagrangian is Vi([(j)*(x), (f>(x)]±) and we represent the 
generating functional for the interacting theory as 



W[J, J*] = exp[i f faVti- E SJ{a f Sj{a) J W free [J, J*}. (22) 

We give just one example of the case where the interaction term in the Lagrangian 
has the nested commutator form, 

Ct = g\[[M*)M*)UM*)U = 9 E ^{x^fixWix) (23) 

all different 

where fa are three different neutral parabose fields of order p — 3. In this case 
we need free Lagrangian terms for each of the three fields, so the free generating 
functional is 



W free [J i ,Jf]=W free [0,0]exp[-i E E j ' d 4 xd 4 yjt\x)A F (x-y)J^(y)}. (24) 

i=j=l ct=l 

The interacting generating functional is 

W[ Ji, J,} = exp[^ J d 3 X E T^JL TAy) \ W free[Ji, ■£]■ (25) 



all different 

Note that for this type of interaction the degree of the fields in the interaction 
Hamiltonian must match the order of the parabose field so that the Green indices 
are "saturated." In addition, the interacting Hamiltonian should commute with the 
field. For even degrees the interaction Hamiltonian will anticommute as pointed out 
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in [13]. This case should not be allowed because it leads to alternating signs of the 
contribution of the interaction energy to widely separated states [15]. 

Now the vacuum matrix elements of the time ordered products of the fields 
4>*{x) and 4>{x) are given by the operations Y/l 8/8 J^ (x) and J2i 5/5J^*(x), re- 
spectively, acting on the generating functional. The anticommutativity of the (p^ 's 
and J*") 's for different values of the Green indices leads to the anticommutativity 
of the quantized fields for different values of the Green index in the vacuum matrix 
elements of the time ordered products, and thus leads to the time ordered products 
of the parabose fields. 

The functional differentiation with respect to J^ (x) can be carried out using 
the following commutation relation, which is analogous to Eq.(7) for the canonical 
formalism, 

jj^J^iv) = ±(2<W - l)jW( y)ll *^ + 8 a , 8(x - y). (26) 

This relation is for a neutral field and the ± refers to parabose (parafermi). For 
a charged field, just omit the contractions that do not connect the field with its 
adjoint. 

As specific illustrations, we evaluate the two-point and four-point Green's func- 
tion for a non-interacting parabose system. We choose to exhibit the results for a 
neutral field, which has all terms that can occur in the case of a charged field. For a 
specific distribution of a charged field and its adjoint use the result for the neutral 
field, but keep only those propagators that connect the charged field and its adjoint. 
Discard those propagators that would have connected the field with the field or the 
adjoint with the adjoint. 



< 01^(^)10 > = EJ-) SJ ^ X2)SJ ^ {xi) w ^--o 



ipA F (x 1 -x 2 ), (27) 



and 

<O|T0(x 4 )0(a;3)0(x 2 )0(x 1 )|O>= 
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J^A" SJ^(x 4 )SJ^(x 3 )SJ^(x 2 )SJ^( Xl ) X |j -° 

= p 2 iA F (xi - x 2 )iA F (x 3 - x 4 ) + p(2 - p)iA F (x 2 - x 4 )iA F (x 1 - x 3 ) 

+p 2 iA F (xi — X4)iA F (x2 — x 3 ). (28) 

In particular, for a charged field, the result is 
< 0\T(f)*(x 4 )(f)*(x 3 )(})(x 2 )(f)(x 1 )\0 >= 

J^W SJ( s )(x A )SJ^(x 3 )SJ^>(x 2 )6J( a >(x 1 ) X Wlj ^=° 
= p(2 — p) iA F (xi — x 3 )iA F (x2 — X4) + p 2 iA F (xi — X4)iA F (x 2 — x 3 ). (29) 

The validity of the above path integral based results can be ascertained through 
the direct evaluation of the time ordered products. In momentum space, the para- 
bose fields are 

1 f d 3 k 



[x 



x 



{2uf' 2 J V^~ k 
d 3 k 



(b k e- ik - x + d\e lk - x ) = ^ + \x) + 0(~)(x), (30) 



(d k e~ ik - x + b\e ik - x ) = (f>< + \x) + (f)* { -\x), (31) 



{2-Kf 2 J V2^~ k 
and these satisfy the vacuum conditions 

0(+)(x)|O>=0* (+) |O>=O. (32) 

Substitution of Eqs.(30, 31) in the time ordered products, subsequent expansion 
of the creation and annihilation operators in terms of Green components, and the 
application of vacuum condition Eq.(32) again lead to the expressions Eq.(26) and 
Eq.(28) for the two and four point Green's functions. 

3.2 Parafermi case 

Next we consider the path integral for parafermions. For a parafermi field we in- 
terchange the roles of commutators and anticommutators relative to the case of 
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parabose fields. Thus the Lagrangian density for a free parafermi field is 



= y2(ip (x) (iy — m)ip^ a '{x)) + c — number. 

a=l 

The generating functional for a free charged parafermi field is 



(33) 



[ f[ Vij (a) *V^ a) exp[-i^ f d 4 x(i> {a) (x)(if - m)ij (a) 



x 



a=l 1 

+f] {a) (x)ij {a) (x) + ^ a \x)7] {a \x))}. (34) 

Now we replace the path integral over commuting fields by the Berezin path inte- 
gral [16] over Grassmann fields and require the Grassmann fields ^ a \x) and ip^ a \x) 
and the external sources ff a ' (x) and rf a ' (x) to anticommute for the same and com- 
mute for different values of the Green indices. The generating functional for a free 
parafermi field is 



WfreeiV, V] = Wf ree [0, 0] exp 



: J2[d A xd 4 y fj {a) (x)S F (x - y)rf a \y) 

a=l 



(35) 



where Sf(x) = (2ir) 4 / d 4 p exp(—ip ■ x){$ + m)(p 2 — m 2 + ie) 1 . In the presence 
of an interaction V(iJj, ip), we modify the generating functional to 



PF [77,77] = exp 



i / d 4 xV{-i 



S 
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Sr](x) ' 5i](x) 



W free [ri,r]]. 



(36) 



These arguments can be extended further to systems in which parafermi and 
parabose fields have mutual interaction. As a specific example of such systems, we 
consider the Yukawa coupling between a charged parafermi and a neutral parabose 
fields of order 3, respectively [13]. These fields obey relative para-Bose commutation 
rules. The interaction term is 



l rr T 



£i = g-[[^(x) : i/j(x)] + - < [i/j(x) : ip(x)]+ > ,<j>(x)\. 

= g Y, tp {a) {x)ilj (f3) {x)<t) M {x). 



(37) 



all different 



n 



The generating functional for the noninteracting system is now 



W free [r),r],J] 



W /ree [0,0,0]exp 



3 

-iJ2 d ^ xdA y n {a \x)S F {x - y)r] (a \y) 

ct=lV 



x exp 



3 f 

: X] / d 4 xd 4 yJ M (x)A F (x - y)J {l \y) 

7=1 



(3* 



which leads to the following expression for the complete generating functional 



W[r),r), J] = exp 



Jd 4 x J2 



s 3 



,13,1 
all different 



ri(<*)(x)fjW(x)jM(x) 



W free [r],fj,J}. (39) 



The propagators for the parafermi system can be obtained by employing the 
commutation relation Eq.(26) with the lower sign. Note that proceeding as in 
the parabose case, we get, for example, the following expression for a four-point 
parafermi propagator 

< 0\T^j{x a )^{x 3 )^{x 2 )^{x 1 )\0 >= 



E 



5 A 



Sf] ( - i4 \x i )5f] {ki) (x 3 )Sri ( - i2) (x 2 )5ri ( - il \x 1 ) 



W 



r)=r)=0 



= p 2 iSp(x 3 — x 2 )i5 , ^(x 4 — xi) — p(2 — p)iSp(x 3 — x\) iSp^x^ — x 2 )- (40) 

One can construct all other generating functionals for para fields, such as the gener- 
ating functionals for connected and for one-particle-irreducible Greens functions, in 
parallel with usual cases, just as we have done above for the time-ordered products. 



4 Remark about the partition function 



The single-particle propagator is related to the partition function [17]. Our result 
Eq.(28) for the two-particle propagator gives a result for the partition function that 
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is the square of the partition function in the special case p = 2 as we expect from 
the result of P. Suranyi [18] that a gas of p = 2 parabosons (or parafermions) is 
equivalent to a mixture of two independent gases of bosons (or fermions). We thank 
A. Polychronakos [19] for pointing out that the different four-point matrix element 
that we gave in an earlier version of this work 

< O|T0*(x 4 )0(a;3)0*(x 2 )0(x 1 )|O >= 

/1\ 4 5 A 

J^W 5J 5 \x A )5J^\x,)5J^\x 2 )5J^\x 1 ) X Wlj= ° 

= p 2 [iA F (xi — x 2 )iA F (x 3 — X4) + iA F (xi — X4)iA F (x 2 — x 3 )]. (41) 

does not have this property. 

5 Remark about normalization of para operators 

We have followed the normalization used by Green [9] when he introduced paras- 
tatistics in which the parameter p that labels the order of the parastatistics does not 
appear in the commutation relations Eq.(3), but p does appear in the normalization 
of the single-particle state Eq.(5). Green's motivation was to keep the usual com- 
mutation relations between the number operator and the annihilation and creation 
operators Eq.(l). One can normalize the single-particle state to one, but then the 
commutation relation becomes 

[[bhb l } ± ,tf m ]_=2p5 lm b{ (42) 

and the number operator commutation relations also acquire the factor p. 

6 Comparison with earlier work 

We are aware of three earlier treatments of parastatistics using path integrals. Y. 
Ohnuki and S.Kamefuchi [20] construct the path integral for parafermi fields using 
paragrassmann variables following the fundamental definition of a path integral as 
the limit of a product of time evolution operators for small time intervals. They 
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do not evaluate this limit and thus do not find an explicit formula for the path 
integral. They also do not give the generating functions that are most useful for 
calculations. Like our construction, their analysis relies on the Green ansatz and 
does not give an intrinsic construction of the path integral. M. Chaichian and A. 
Demichev [21] and A. Polychronakos [22] both give an intrinsic definition of the 
path integral that does not rely on the Green ansatz. Their construction of the 
path integral is a first quantized analog of using the Green trilinear commutation 
relations without making use of the Green ansatz. They define the propagator for 
paraparticles using the propagator for distinguishable particles together with the 
matrix elements of the permutation operator. They point out that these matrix 
elements are "weights for each topological sector of the path integral" that are 
simply positive or negative integers or zero. The result of Polychronakos that, in his 
notation, S p F(P) = (—l) p S p B(P), follows in our case from the ± sign in Eq.(7,26). 

7 Summary and conclusions 

We have shown that parastatistics fields can be quantized using path integrals. All 
the results are, as we expect, the same as found with the canonical formalism. Up to 
now we have not found a computational advantage in using either of these formalism 
for parastatistics. We explicitly calculated the two- and four-point functions for free 
parabose and free parafermi theories and also gave the generating functionals for 
both free and interacting theories. It would be interesting to formulate the path 
integral for parastatistics fields using only a definition of the path integral that 
uses the trilinear commutation relations for parafields without reference to Green's 
ansatz. 
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